function F1= delayFFT(F,T,fs)

%M debe ser par
M=length(F);
delay=fs*T;
f= [0:(M/2-1),-M/2]'; %Vector columna
aux=-2*pi*f*mod(delay/M,1);
Wb=cos(aux)+1i*sin(aux);
Wb(M/2+1)=cos(aux(end));
F1=F(1:(M/2+1)).*Wb;
F1=[F1;conj(F1((M/2):-1:2))];